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Abstract:  In  the  development  of  the  guidance  and 

control  packages  for  unmanned  vehicles,  it  is  highly 
desirable  to  have  inertial  measurement  sensors  which 
are  small,  inexpensive,  low  power,  reliable  and 
accurate.  New  technological  advances  in  the  design 
and  construction  of  micro  inertial  sensors,  such  as 
accelerometers  and  gyroscopes,  have  much  promise  in 
providing  small,  inexpensive,  and  low  power  devices; 
however,  much  improvement  in  the  reliability  and, 
especially,  the  accuracy  of  these  micro  devices  is  still 
necessary.  Further  major  improvements  in  these  two 
properties  will  probably  not  be  accomplished  in  the  near 
future,  thus  it  will  be  necessary  to  use  special  signal 
processing  methods  to  provide  the  accuracy.  One  way 
which  has  been  proposed  to  improve  the  accuracy,  and 
concurrently  the  reliability,  of  micro  sensors  is  to  use 
many,  perhaps  one  hundred  or  more,  micro  sensors  on  a 
single  chip  (or  a few  chips)  and  using  statistical 
methods  to  combine  the  individual  outputs  of  these 
sensors  to  provide  an  accurate  measurement.  One 
method  of  performing  such  a combination  is  through  an 
extended  Kalman  filter  (EKF).  A standard  application 
of  an  EKF  to  an  array  of  gyroscopes  would  involve  at 
least  six  state  equations  per  gyroscope  and  the  number 
of  covariance  equations  would  be  in  the  order  of  the 
square  of  the  product  of  six  times  the  number  of 
gyroscopes.  Obviously,  the  ’curse  of  dimensionality’ 
very  quickly  limits  the  number  of  sensors  (gyroscopes) 
which  can  be  used.  Even  if  the  EKF  for  each  individual 
gyroscope  is  uncoupled  from  the  rest,  the  number  of 
covariance  equations  is  of  the  order  of  the  number  of 
gyroscopes  times  six  squared.  This  can  still  lead  to  a 
formidable  computational  burden.  In  this  paper,  a new 
technique  of  applying  an  EKF  to  this  problem  of 
combining  many  sensors  is  proposed.  By  using  the 
common  nominal  model  for  each  of  the  micro  sensors 
and  developing  a single  EKF,  improved  accuracy  is 
achieved  by  a single  EKF  with  the  dimension  of  one 
sensor.  For  cases  in  which  the  micro  sensors  are 
corrupted  by  correlated  noise  (between  the  sensors)  an 
artificial  neural  network  could  be  added  to  the  EKF 
dynamics  to  track  the  noise.  Simulated  examples  will 
be  discussed. 

Summary7:  Recent  advances  in  the  technology  of 

microelectromechanical  systems  (MEMS)  has  led  to 
optimistic  predictions  of  the  potential  applications  of 
various  micro  devices.  One  family  of  devices  which 
appears  to  have  wide  applicability  in  many  areas  is  the 
family  of  inertial  sensors,  including  gyroscopes  and 


accelerometers.  Micro- sensors,  because  of  their  small 
size,  often  have  only  moderate  accuracy  as  compared 
to  the  accuracy  of  full-size  sensors.  Applications  which 
demand  sensors  with  small  size  and  low  cost  as  well  as 
high  accuracy,  will  require  signal  processing  methods. 
Because  of  the  expected  low  cost  of  the  micro-sensors, 
one  approach  to  increased  accuracy  is  to  use  a large 
number  of  micro-sensors  to  measure  the  same  quantity 
and  then  use  statistical  methods  to  combine  the  many 
low  accuracy  measurements  to  generate  a single  high 
accuracy  measurement.  For  many  years  the  extended 
Kalman  filter  (EKF)  has  been  used  to  process  output 
signals  from  inertial  sensors  to  produce  more  accuracy 
than  is  possible  from  the  raw  output  signals,  thus  it  is 
natural  to  assume  that  the  EKF  can  be  readily  adapted 
to  the  many  micro-sensor  problem.  In  this  paper,  a 
technique  is  proposed  in  which  an  arithmetic  average  of 
the  many  outputs  from  a set  of  micro-gyroscopes  is 
input  to  a single  average  EKF  Under  certain  simple 
assumptions,  this  can  be  shown  to  produce  the  desired 
result.  Using  a model  of  a real  micro-gyroscope  this 
technique  was  simulated  for  a set  of  ten  micro- 
gyroscopes  and  the  preliminary  results  indicate  that  this 
method  has  promise.  Continuing  exploration  of  this 
method  will  include  adaptive  EKF  methods  using  a 
neural  network  based  EKF  which  will  allow  the  filter  to 
track  correlated  errors. 

Introduction.  Recent  developments  in  micro- 
electromechanical systems  (MEMS)  technology  have 
led  to  optimistic  predictions  for  the  use  of  micro 
devices  in  a wide  range  of  applications.  It  appears  that 
micro-sensors,  in  particular,  have  immediate 
applicability  in  many  fields  from  medical  implants  to 
automobiles  to  aerospace  vehicles.  Small  size  plus  their 
probability  of  becoming  very  inexpensive  make  their 
potential  range  of  applications  almost  limitless. 
However,  in  spite  of  small  size  being  a major  advantage 
of  micro-sensors,  it  can  also  be  a major  disadvantage. 
In  particular,  in  many  cases,  the  accuracy  of  a sensor 
can  be  controlled  much  easier  if  the  sensor  has  large 
dimensions,  because  accuracy  is  often  a function  of 
dimensional  ratios  and  controlling  such  ratios  by 
machining  is  easier  when  dimensions  are  larger.  On  the 
other  hand,  if  the  micro-sensors  can  be  made  much 
more  inexpensively  than  larger  machined  sensors  then 
this  disadvantage  may  be  overcome  by  using  many 
relatively  inaccurate,  but  also  inexpensive,  sensors  as 
opposed  to  one,  or  a few,  accurate,  but  much  more 
expensive,  sensors.  Note  that  the  use  of  many  micro- 
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sensors  measuring  the  same  quantity  can  also  provide 
reliability  through  redundancy  at  a reasonable  cost 
This  latter  point,  while  of  considerable  interest,  will  not 
be  considered  in  this  paper.  The  problem  of  interest 
here  is  that  of  combining  the  outputs  of  several  micro- 
sensors, all  measuring  the  same  quantity,  so  that  the 
accuracy  of  the  combination  exceeds  the  accuracy  of 
the  individual  micro-sensors  but  without  requiring  an 
extraordinary  amount  of  computation.  This  is  what  will 
be  referred  to  as  the  signal  processing  problem  for 
micro  inertial  sensors. 

The  signal  processing  problem:  A sensor  is  a device 
that  measures  a physical  quantity  and  produces  a 
corresponding  output,  typically  an  electrical  quantity, 
which  is  related  in  a known  way  to  the  physical 
quantity.  In  practice,  the  measurement  of  the  physical 
quantity  is  corrupted  by  noise  and  the  relationship 
between  the  physical  quantity  and  the  corresponding 
output  is  corrupted  by  another  noise.  The  accuracy  of 
the  sensor  is  dependent  upon  the  magnitudes  of  these 
noises.  The  accuracy  of  a sensor  might  be  improved  by 
better  understanding  of  the  physics  by  which  the  sensor 
measures  the  physical  quantity  and  how  it  produces  the 
corresponding  output  and  improving  the  process  by 
which  the  sensor  is  manufactured  or  by  applying 
appropriate  signal  processing  techniques  to  the  output 
signal,  In  the  area  of  micro-sensors,  significant 
accuracy  improvements  in  the  manufacturing  process 
will  require  increased  costs  and/or  larger  geometries, 
both  of  which  will  tend  to  neutralize  the  advantages  of 
such  devices.  Thus  signal  processing  techniques 
applied  to  multiple  sensor  configurations  are  being 
examined  and  are  expected  to  improve  the  accuracy  of 
existing  micro-sensors.  This  is  not  to  imply  that 
improved  manufacturing  processes  are  not  being 
developed;  however,  currently,  appropriate  processing 
of  micro-sensor  outputs  appears  to  be  the  quickest  and 
most  feasible  method  of  improving  sensor  accuracy. 

The  simplest  signal  processing  concept  to  improve 
micro-sensor  accuracy  is  to  manufacture  one  (or  a few) 
chips  with  a total  of  many  micro-sensors  all  of  which 
measure  the  same  physical  quantity.  The  outputs  of  all 
of  the  micro-sensors  are  then  simply  averaged 
arithmetically.  If  the  noises  associated  with  the  outputs 
of  the  several  micro-sensors  are  additive,  mutually 
independent  and  zero  mean,  then  the  standard  deviation 
of  the  error  of  the  arithmetic  average  is  reduced  by  the 
square  root  of  the  number  of  output  signals.  Thus,  if 
one  hundred  micro-sensors  are  used,  it  is  equivalent  to 
replacing  them  with  a single  micro-sensor  with  an  error 
with  one-tenth  the  standard  deviation.  However,  if  the 
errors  in  the  output  signals  are  non-zero  mean  (across 
the  ensemble  of  micro-sensors)  or  if  the  errors  are 
correlated  across  the  micro-sensors,  then  arithmetic 
averaging  may  not  be  effective.  That  is,  the  averaging 
will  reduce  the  uncorrelated  component  of  the  noise  but 


may  have  little  effect  on  the  correlated  component  of 
the  noise.  A correlated  component  of  the  noise  in  the 
output  signals  is  to  be  expected  if  many  of  the  sensors 
are  on  a single  chip.  For  example,  any  alignment  errors 
in  the  chip  manufacturing  process  may  generate 
measurement  errors  in  one  micro-sensor  which  are 
strongly  correlated  with  all  of  the  measurement  errors 
on  the  other  micro-sensors.  If  the  manufacturing 
process  generates  alignment  errors  which  are  correlated 
from  chip  to  chip,  for  example  all  the  chips  come  from  a 
single  wafer,  then  this  correlation  may  extend  to  all  of 
the  micro-sensors.  In  this  case,  simple  arithmetic 
averaging  may  not  be  effective  as  a signal  processing 
technique. 

In  the  application  of  a standard  sized  inertial  sensor,  it  is 
quite  common  to  process  the  output  of  the  sensor 
through  a Kalman  filter,  or  more  likely  through  an 
extended  Kalman  filter  (EKF),  to  improve  the  accuracy 
of  the  measurement.  If  a number  of  inertial  sensors,  for 
example,  three  gyroscopes  and  three  accelerometers, 
are  combined  on  a single  platform  it  is  common  that  a 
single  extended  Kalman  filter  is  used  to  generate 
estimates  of  the  state  of  the  platform  on  which  the 
sensors  are  mounted.  These  estimates  are  generally 
much  improved  over  the  measurements  which  are  taken 
directly  from  the  sensors.  A similar  concept  could  be 
used  to  generate  an  estimate  of  the  common  quantity 
being  measured  by  the  many  micro-sensors.  In 
particular,  micro  inertial  sensors  have  well  developed 
mathematical  models  and  are  well  suited  to  the 
application  of  an  extended  Kalman  filter.  The  difficulty 
with  applying  this  concept  to  a set  of  many  micro 
inertial  sensors  is  that  the  dimensionality  of  the  EKF 
may  become  excessive.  For  example,  if  a single 
gyroscope  is  modeled  by  a set  of  six  coupled  dynamic 
equations,  then  six  hundred  equations  would  be  needed 
to  model  one  hundred  gyroscopes.  The  resultant  EKF 
which  uses  this  dynamic  model  would  have  six  hundred 
states  and  its  covariance  equation  would  have  between 
two  thousand  and  one  hundred  eighty  thousand 
(180,000)  equations  depending  upon  the  coupling 
between  the  various  gyroscopes.  Obviously,  this 
method  of  combining  measurements  is  limited  by  the 
number  of  micro  gyroscopes  which  are  being  used. 

The  method  proposed  in  this  paper  for  combining  the 
measurements  from  many  micro-sensors  is  to  first 
generate  the  arithmetic  average  of  the  outputs  of  the 
sensors  and  then  process  this  average  output  through  an 
average  EKF.  Consider  the  following  standard  Kalman 
filter  problem  of  estimating  the  state  of  a system  defined 
by  the  vector-matrix  equation 
ch 

— = Ax  + Bu  + w ( 1 ) 

dl 

where  x is  the  system  state  vector,  u is  the  system 
input  vector,  and  w is  a zero  mean  white  noise  vector 
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with  covariance  Q . A and  B are  the  matrices  which 
define  the  system.  The  state  vector  is  measured 
through  a noisy  linear  transformation 

y rr  CX  + V (2) 

where  C defines  the  transformation  and  v is  a zero 
mean  white  noise  with  covariance  matrix  R.  Suppose 
that  the  system  state  can  also  be  measured  from  a set  of 
Another  systems  defined  by 
dx : 

= Ax,-i-Bii  + w , / = 1,2,  — , 

dt  1 

where  x is  the  state  vector  and  w is  the  white  noise 
/ / 

vector  of  the  /th  system.  The  states  of  these  systems 

are  measured  through  a set  of  N measurement 

equations 

y /-=Cx/-+v/-,  / = 1,2,  — , N 


where  v,  is  the  white  noise  of  the  /th  measurement. 
An  arithmetic  average  of  the  outputs  is  formed  as 


- 1 r 

y=v;s, 


( 


h = C| 


1 


N 


— Z X. 

Ni  = 1 ' 


N 


N 

I v. 
/ = 1 ' 


or 


y = Cx  + v 


(3) 


Now  averaging  the  N differential  equations  results  in 
the  average  differential  equation 
dx 

— = Ax  + Bu  4-  w (4) 

dt 

In  the  stochastic  sense  equation  (3)  is  identical  to 
equation  (2)  and  equation  (4)  is  identical  to  equation 
(1),  except  that  the  noise  w has  a covariance  of 
Q / N and  the  noise  v has  a covariance  of  R IN. 


Thus  if  estimating  the  system  response  to  the  input  u is 
of  primary  interest,  the  effect  of  the  white  noises  is 
reduced  if  the  average  system  output  is  used  to  drive 
the  Kalman  filter  rather  than  the  measurement  in 
equation  (2).  In  the  following  example,  the  problem  of 
estimating  the  output  of  a set  of  ten  micro-gyroscopes 
all  measuring  the  same  angular  rate  using  a single  EKF 
was  simulated  and  compared  with  the  method  of 
applying  ten  EKFs  with  output  averaging 


Mathematical  model  of  a micro-gyroscope:  In  this 
section,  the  mathematical  model  of  a micro  gyroscope  is 
presented.  The  parameters  used  later  in  the  simulation 
are  taken  from  an  actual  experimental  micro-gyroscope 
developed  by  Irvine  Sensors,  thus  the  simulation 
represents  the  real  world  as  nearly  as  we  were  able  to 
model  it.  The  gyroscope  is  composed  of  a small  rigid 
plate  attached  to  a housing  through  a set  of  four  thin 
orthogonal  support  wires.  A set  of  orthogonal  axes  is 
fixed  to  the  plate  with  the  origin  at  the  center  of  mass  of 
the  plate  and  the  three  axes  are  aligned  with  the 
principal  axes  of  the  plate.  The  x-axis  and  the  y-axis  lie 
in  the  plane  of  the  plate  and  are  assumed  to  be  aligned 
with  the  support  wires.  In  an  application,  the  housing  is 
attached  firmly  to  a much  larger  rigid  body  (perhaps  an 


automobile),  thus  the  plate  can  rotate,  to  a limited 
extent,  about  all  three  axes  relative  to  the  housing  and, 
thus  also  rotates  relative  to  the  larger  rigid  body. 
However,  it  is  assumed  that  that  any  translation  relative 
to  the  larger  rigid  body  is  negligible  Because  the  axes 
rotate  with  the  plate,  the  products  of  inertia  of  the  plate 
in  this  coordinate  system  are  zero.  The  plate  is  forced 
to  oscillate  at  a constant  amplitude  and  frequency  about 
the  2-axis  by  a sinusoidal  torque  applied  by  an  electric 
field.  The  x-axis  is  the  input  axis  through  which  the 
input  angular  velocity  of  the  large  rigid  body  is  coupled 
into  the  plate  through  the  support  wires.  The  y-axis  is 
the  output  axis  of  the  gyroscope  from  which  a 
measurement  of  the  input  angular  velocity  is  obtained  as 
some  function  of  the  periodic  motion  induced  about  the 
>’-axis  by  the  angular  motions  about  the  x-axis  and  the 
2-axis.  Under  these  conditions,  Euler’s  equations  of 
motion  can  be  used  to  develop  a dynamic  model  of  the 
plate.  The  plate  is  symmetric  in  the  x-y  plane  so  that 
the  moments  of  inertia  about  the  x-axis  and  the  y- axis 
are  equal.  The  moment  of  inertia  about  the  2- axis  is 
much  larger  than  these  two.  We  define  the  angular 
position  of  the  plate  from  its  ‘rest’  position  by  the 
angles  6 x ? 0 v , and  6 _ measured  about  the  x,  y,  and  2 

axes.  A sinusoidal  torque,  T(1),  is  applied  to  the  plate 
about  the  2-axis  producing  a periodic  motion  about  the 
2-axis.  The  other  torques  about  the  2— axis  are  a 
damping  torque  and  a spring  torque,  both  due  to  the 
mechanical  properties  of  the  supporting  wires  along  the 
x-axis  and  the  j'-axis.  The  applied  torques  about  the  y- 
axis  and  the  x-axis  are  a damping  torque  and  a spring 
torque  due  to  the  supporting  wires  and  the  damping 
coefficients  and  spring  constants  are  assumed  to  be 
equal  about  each  axis  because  of  the  symmetry  of  the 
plate  in  the  x-y  plane.  When  the  housing  is  rotated  at 

an  angular  rate  of  Qx  rad/sec  about  the  x-axis,  a 
torque  Tx(t)  is  transmitted  to  the  plate  by  the  support 

wires  and  consists  of  a damping  torque  and  a spring 
torque 

t 

Cx  (Qx  -0J,  and  Kx  (\  Q(l  ')dt  '-0J 

0 

where  C and  Kx  are  the  damping  coefficient  and  the 

spring  constant  about  the  x-axis.  The  input  rate  Q v 
induces  a motion  in  the  y- axis  and  Q v is  determined 
from  the  measurement  of  the  angular  motion  about  the 
>?-axis.  Based  on  these  conditions  , Euler’s  equations 
can  be  written  as 

t 

0,  + ot0x  + aa6x  + A0r0.  =a,  Qv  + aAo.Jt')dt' 

0 

0t  + ax0r  + ao0v-A9x0:  = 0 
0.  + bt  0.  + b09.  = Tu  sir \a)J 
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where  the  coefficients  are  determined  from  the  various 
physical  parameters  of  the  gyroscope.  The  desired 
output  of  the  gyroscope  is  the  input  angular  rate  output 
Qv(/) . This  is  obtained  by  measuring  and  processing 

the  angular  rate  9 (/) . 

Steady-State  Operation:  Assuming  that  the  rate  of 

change  of  the  angular  rate  input  Qv(/)to  the 
gyroscope  is  slow  relative  to  the  rates  of  change 
induced  by  the  applied  torque  70  sin#;0/  and  that  the 
nonlinear  effects  in  the  x and  y equations  are  relatively 
small  phasor  analysis  can  be  used  to  determine  a 

steady-state  value  of  the  output  6 (Oas 

y 

i _ 

0 (t)-A'IG  (co)G  (ft)  )ft)~Q  sin  ft)  / 
y 1 0 rv  (T  yy  CT  0 x 0 

where 

i , ! , 

V'W  4 'Vi* 

and 


V'V3 


Note  that  the  phase  shift  of  the  phasor  is  not  needed 
and  has  been  ignored.  Note  that  this  equation  could  be 
written 


9 (/)  - KQ  (z)sinft;  1 
y x 0 


where  K is  a constant.  Thus  6 (!)  is  sine  wave 

3’ 


sin  ftW  modulated  by  the  applied  angular  rate  Q (/) . 

The  angular  rate  can  thus  be  obtained  by  a simple 
amplitude  modulation  (AM)  demodulator. 


State  Equations  for  the  Gyroscope  For  a Kalman 
filter  the  equations  for  the  gyroscope  must  be  written  in 
the  form  of  state  equations.  To  put  the  three  equations 
into  a state  equation  form,  let 

*2  = 0X 

x3  = 6 y -v4  = .v,  = 0V 

*5  = = j5  = 9 z 

then  the  state  equation  representation  for  the  gyroscope 
is 

i,  = ~a{]x2  +a0  Qv 

x2  = a",  - ci\X2  - Av4 x0  +a ,QV 

i,  = a-4 

i4  =-a{) a,  -o,jr4  + Ax :2x(, 
x5  = a-6 

X(,  =-b0x5-blx(>  + r0sm(o0t 


where  the  variable  A',  has  no  specific  physical  meaning. 

Two  additional  equations  can  provide  a simple 
demodulation  of  the  variable 


■r4(t)  = 0yU) 


to  generate  Q 

x 


(/) . These  equations  are 


A' 7 — A"s 


i8  - -6*0a7  -l\Xz+c0xa  sinftV 
where  the  coefficients  were  chosen  to  provide  low  pass 
filtering  with  a corner  frequency  at  CO{)  / 10  rad/sec. 
The  output  measurement  is  now  the  modulator  output, 
that  is, 

-ft)  = X7(t) 

Ideally,  in  steady  state,  this  is  related  to  Q (/)  by  the 

relationship  developed  in  the  steady-state  section  and 
the  demodulator  gain  by 

x 2 7 

Design  Coefficients  The  design  coefficients  for  the 
gyroscope  simulated  in  this  paper  are 

a,  =1.8621  xlO4 
a0  =3.4483  x I08 
/>,  = 9.3750x10 
h0=  3.3750  xlOx 
A = 4.5172 

and 

co0  = 2/r  x 3000  = 1.8850  x 1 04  rad  / sec 
which  is  near  the  resonant  frequency  of  the  linear  part 
of  the  y-axis  equation.  The  demodulator  coefficients 
are 

c j = 1200/t  and  c = (600^)"' 

The  magnitude  /()  of  the  forcing  function  was  chosen 
to  be 

T0  = 1.5  x 108 

which  generates  a maximum  steady  state  motion  about 
the  x-axis  of  a fraction  of  a radian  The  resulting  state 
equations  using  these  coefficients  was  used  as  the  Truth 
model  for  the  simulations  which  are  discussed  later. 


Kalman  Filter  Equations  To  generate  a set  of 
Kalman  filter  equations,  the  following  modification  to 
the  gyroscope  state  equations  had  to  be  made.  Let: 

Q v = a' 2 + n] 

where  n j is  the  unknown  difference  between  the  input 

fiv  and  the  state  x2  —0X.  The  resulting  state 
equations  can  now  be  rewritten: 
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*2  =xi-Ax4xb+aint 

i,  = A'4 

i4  = -a0x3  - arx4  + Ax2x6  + n4 


*6  = ~hox5  - htx6  + y0  sin  6)J  + n6 

with  the  demodulator  equations 
X7  = x8 

x8  = — c0jc7  - cxx%  + cQxA  sin  co{)1  + w8 

The  measurement  equation  is  given  by 

z = x7  +v 

The  quantities  nx , n4 , , n8  and  V are  random 

quantities  which  in  the  development  of  the  extended 
Kalman  filter  (EKF)  are  assumed  to  be  white  noises. 
As  almost  four  decades  of  usage  has  shown  the  EKF  is 
robust  to  inaccuracies  in  the  assumptions  on  the  system 
noises.  For  convenience,  the  gyroscope  model  is  now 
rewritten  in  the  vector-matrix  form: 

x = f (x)  + Gn  + g(x)  sin  co0t 


- = hrx  + v 
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and 

f(x)  = Fx  + f,(x) 

where 


F = 


and 
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The  noise  statistics,  which  are  normally  determined 
empirically,  are  chosen  to  have  the  following  properties, 
R is  a nonzero,  positive  scalar  and  Q is  a nonzero, 
positive  definite  4x4  matrix  assumed  to  be  diagonal 
and  given  by: 


where 
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The  extended  Kalman  filter  equations  can  now  be 
written  as 


x = f(x)-hk(r-h7x)  + g(x)*sin^0/ 

where  X , the  output  of  the  Kalman  filter,  is  the 
estimate  of  the  combined  state  of  the  gyroscope  and  the 
demodulator.  The  estimate  of  the  applied  torque 
Qv(/)  is  given  by  the  state  estimate  J?7 (7 ) . 

k - Ph  / R is  the  Kalman  gain  vector  and  the  error 
covariance  P is  calculated  by  the  equation 
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_ Phh7P 

P = FP  + PF + GQG' 

R 

where 

F = F + F, 

and 
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0 Ax,  0 0 0 Ax,  0 0 

F,  = ( 
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00000000 
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These  equations  were  used  to  simulate  the  EKFs  used 
in  the  simulations. 

Simulation  results:  The  technique  described  in  this 

paper  of  averaging  the  output  measurements  from 
several  micro-gyroscopes  and  processing  the  average 
output  in  a single  EKF  to  estimate  the  common  angular 
rate  applied  to  each  micro-gyroscope  was  simulated 
using  ten  micro-gyroscopes  with  the  micro-gyroscope 
model  described  earlier  The  applied  angular  rate  input 
Q_v(/)  was  a sine  wave  with  a maximum  amplitude  of 
one-tenth  rad/sec  and  a frequency  of  one  hundred  Hertz 
.A  measurement  noise  with  a standard  deviation  of 
twenty  per  cent  of  the  maximum  of  the  steady-state 
value  (,012)  of  the  output  signal  of  the  truth  model 
driven  by  this  input  was  used  to  cormpt  the  output  of 
each  micro-gyroscope.  The  output  of  this  simulation 
was  compared  to  the  simulated  output  of  the  truth 
model.  The  results  are  shown  in  the  figure  below  in 
which  the  solid  line  is  the  noise  free  output  and  the 
dotted  line  is  the  estimate.  The  standard  deviation  of 
the  time  average  of  the  steady  state  error  is 

(23672)10  4 . This  is  an  improvement  of  more  than 
ten  times  over  the  given  measurement  noise.  With 

simple  averaging,  an  improvement  of  -y/TO  would  be 
expected. 


Estimate  of  Output  vs  Noise-free  Output  (1  EKF  Case) 


Time  m Seconds 


The  technique  of  developing  individual  estimates  of  the 
applied  angular  rate  using  one  EKF  for  each  of  the  ten 
micro-gyroscopes  and  averaging  the  outputs  was  also 
simulated  and  compared  to  the  truth  model.  These 
results  are  shown  in  the  figure  below  in  which  the 
dotted  line  is  the  noise-free  output  and  the  solid  line  is 
the  estimate  The  standard  deviation  of  the  time 

average  of  the  steady  state  error  is  (2.1 058)1 0 4 . 

Estimate  of  Output  vs  Noise-free  Output  (10  EKF  Case) 
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Finally,  the  two  multiple  sensor  techniques  were 
compared  to  each  other.  The  error,  compared  to  the 
truth  model,  was  generated  for  each  and  the  difference 
between  these  two  errors  is  shown  in  the  figure  below. 


Time  in  Seconds 


It  is  apparent  that  the  two  different  signal  processing 
methods  are  very  comparable.  There  seems  to  be  only  a 
slight  degradation  of  the  estimate  when  one  EKF  is 
used  as  opposed  to  ten  EKFs.  On  the  other  hand,  the 
amount  of  computation  is  considerably  reduced. 

Conclusions:  In  this  paper  a method  of  improving  the 
accuracy  of  micro  inertial  sensors  is  proposed  This 
technique  uses  multiple  sensors  with  advanced  signal 
processing  techniques  to  combine  the  many  signals  into 
a single  output.  Two  such  techniques  are  discussed,  the 
first  is  very  computationally  intensive  and  the  second  is 
moderately  computationally  intensive.  In  order  to  be 
useful,  the  outputs  generated  by  these  techniques 
should  be  significantly  more  accurate  than  the  output  of 
a single  sensor  or  even  of  the  arithmetic  average  of  the 
outputs  of  many  sensors.  Some  preliminary  simulation 
results  indicate  that  both  techniques  improve  the  sensor 
accuracy  significantly;  however,  the  computational 
intensity  of  the  first  technique  increases  exponentially 
with  the  number  of  sensors,  which  precludes  it  from 
being  used  with  a very  large  number  of  sensors.  In  the 
second  technique,  the  computational  burden  is 
independent  of  the  number  of  sensors  and  if  it  has 
comparable  accuracy  to  the  first  technique  then  is  far 
preferable.  Early  comparisons  of  the  two  techniques 
indicates  that  the  accuracies  are  comparable. 

This  research  is  still  in  a preliminary  stage  and  definitive 
conclusions  are  premature;  however  the  early 
simulations  give  results  that  are  promising  enough  to 
encourage  further  work  in  this  area. 


